The angular size of the Cepheid I Car: a comparison of the 
interferometric and surface brightness techniques 
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ABSTRACT 

Recent interferometric observations of the brightest and angularly largest classical Cepheid, 
£ Carinae, with ESO's VLT Interferometer (VLTI) have resolved with high precision the variation 
of its angular diameter with phase. We compare the measured angular diameter curve to the one 
we derive by an application of the Baade-Wesselink type infrared surface brightness technique, 
and find a near-perfect agreement between the two curves. The mean angular diameters of ^ Car 
from the two techniques agree very well within their total error bars (1.5%), as do the derived 
distances (4 %). This result is an indication that the calibration of the surface brightness relations 
used in the distance determination of far away Cepheids is not affected by large biases. 

Subject headings: Stars: distances - Stars: fundamental parameters - Stars: variables: Cepheids - Stars: 
oscillations - Techniques: interferometry - extragalactic distance scale 



1. Introduction 

Cepheid variables are fundamental objects for 
the calibration of the extragalactic distance scale. 
Distances of Cepheids can be derived in at least 
two different ways: by using their observed mean 
magnitudes and periods together with a period- 
luminosity relation, or by applying a Baade- 
Wesselink (hereafter BW) type technique to de- 
termine their distances and mean diameters from 



^European Southern Observatory, CasiUa 19001, Santi- 
ago 19, Chile 

^ Observatoire Midi-Pyrenees, UMR 5572, 14, avenue 
Edouard BeUn, F-31400 Toulouse, France 

^ Astrophysikalisches Institut Potsdam, An der Storn- 
warte 16, D-14482 Potsdam, Germany 

^Universidad de Concepcion, Departamento de Fisica, 
CasiUa 160-C, Concepcion, Chile 

^ Space Telescope Science Institute, 3700 San Martin 
Drive, Baltimore, MD 21218, USA 

•^GEMINI, UMR 6203, Observatoire de la Coto d'Azur, 
Avenue Copernic, F-06130 Grasse, France 

'^LESIA, Observatoire de Paris-Meudon, 5, place Jules 
Janssen, F-92195 Meudon Cedex, Prance 



their observed variations in magnitude, color and 
radial velocity. This latter technique has been 
dramatically improved by the introduction of the 
near-infrared surface brightness method (hereafter 
IRSB) by Welch (1994), and later by Fouque & 
Gieren (1997) who calibrated the relation between 
the y-band surface brightness and near-infrared 
colors of Cepheids. For this purpose, they used 
the observed interferometric angular diameters of 
a number of giants and supergiants bracketting 
the Cepheid color range. This method has been 
applied to a large number of Galactic Cepheid 
variables, for instance by Gieren et al. (1997), 
Gieren et al. (1998), and Storm et al. (2004). 

Applying the surface-brightness relation de- 
rived from stable stars to Cepheids implicitly as- 
sumes that the relation also applies to pulsating 
stars. The validity of this assumption can now 
be addressed by comparing direct interferometric 
measurements of the angular diameter variation 
of a Cepheid to the one derived from the IRSB 
technique. It has recently been shown by Kervella 
et al. (2004), hereafter K04, that the VLT In- 
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terfcromctcr on Paranal is now in a condition to 
not only determine accurate mean angular diam- 
eters of nearby Cepheid variables, but to follow 
their angular diameter variations with high preci- 
sion. Using the Palomar Testbed Interferometer, 
Lane et al. (2000, 2002) resolved the pulsation of 
the Cephcids ( Gem and 77 Aql as early as 2000, 
but the comparison we present in this letter is 
the first where error bars on the derived distance 
and linear diameter arc directly comparable at a 
few percent level between the interferometric and 
IRSB techniques. 

The star that wo will discuss in this letter, 
£ Car, is the brightest Cepheid in the sky. Its long 
period of about 35.5 days implies a large mean 
diameter, which together with its relatively short 
distance makes it an ideal target for resolving its 
angular diameter variations with high accuracy. In 
this paper, wc will compare the interferometrically 
determined angular diameter curve of £ Car with 
that determined from the IRSB technique, and we 
will demonstrate that the two sets of angular di- 
ameters are in excellent agreement. Based on the 
available high-precision angular diameter and ra- 
dial velocity curves for this star, we will also derive 
a revised value of its distance and mean radius. 

Several authors (Sasselov & Karovska 1994; 
Marengo et al. 2003, 2004) have pointed out po- 
tential sources of systematic uncertainties in the 
determination of Cepheid distances using the in- 
terferometric BW method. In particular, imper- 
fections in the numerical modeling of Cepheid at- 
mospheres could lead to biased estimates of the 
limb darkening and projection factor. We will dis- 
cuss the magnitude of these uncertainties in the 
case of £ Car. 

2. Interferometric observations 

The interferometric observations of £ Car were 
obtained with the VLT Interferometer (Glinde- 
mann et al. 2000), using its commissioning in- 
strument VINCI (Kervella et al. 2000, 2003) and 
0.35 m test siderostats. This instrument recom- 
bines the light from two telescopes in the infrared 
K band (2.0- 2.4/im), at an effective wavelength 
of 2.18 /xm. A detailed description of the interfer- 
ometric data recorded on £ Car can be found in 
K04. 

The limb darkening (LD) models used to de- 



rive the photospheric diameters from the fringe 
visibilities were taken from Claret (2000). The 
correction introduced on the uniform disk (UD) 
interferometric measurements by the limb darken- 
ing is small in the K band: for £ Car, we determine 
k = ^ud/^ld = 0.966. Considering the magnitude 
of this correction, a total systematic imccrtainty 
of ±1 % appears reasonable. However, until the 
limb darkening of a sample of Cepheids has been 
measured directly by interferometry, this value re- 
lies exclusively on numerical models of the atmo- 
sphere. This is expected to be achieved in the 
next years using for instance the longest baselines 
of the VLTI (up to 202 m) and the shorter J and 
H infrared bands accessible with the AMBER in- 
strument (Petrov et al. 2000). 

The LD correction is changing slightly over the 
pulsation of the star, due to the change in effec- 
tive temperature, but Marengo et al. (2003) have 
estimated the amplitude of this variation to less 
than 0.3% peak to peak in the H band (for the 
10 days period Cepheid C Gem). It is even lower 
in the K band, and averages out in terms of rms 
dispersion. As a consequence, we have neglected 
this variation in the present study. 

The limb-darkened angular diameter measure- 
ments are listed in Table 1. Two error bars are 
given for each point, corresponding respectively 
to the statistical uncertainty (internal error) and 
to the systematic error introduced by the uncer- 
tainties on the assumed angular diameters of the 
calibrator stars (external error) . The phases given 
in Table 1 are based on the new ephemeris derived 
in Sect. 4. These measurements were obtained dur- 
ing the commissioning of the VLTI, and part of 
them are affected by relatively large uncertain- 
ties (3-5%) due to instrumental problems. How- 
ever, the precision reached by VINCI and the test 
siderostats on this baseline is of the order of 1% 
on the angular diameter, as demonstrated around 
the maximum diameter phase. 

3. The infrared surface brightness tech- 
nique 

The IRSB technique has been presented and 
discussed in detail in Fouque & Gieren (1997) 
(hereafter FG97). In brief, the angular diameter 
curve of a given Cepheid variable is derived from 
its V light and {V — K) color curve, appropri- 
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atcly corrected for extinction. It is then combined 
with its Hnear displacement, which is essentially 
the integral of the radial velocity curve. A linear 
regression of pairs of angular diameters and lin- 
ear displacements, obtained at the same pulsation 
phases, yields both the distance, and the mean ra- 
dius of the star. 

While there are several sources of systematic 
uncertainty in the method, as discussed in Gieren 
et al. (1997), one of its great advantages is its 
strong insensitivity to the adopted reddening cor- 
rections, and to the metallicity of the Cepheid 
(Storm et al. 2004). With excellent observa- 
tional data at hand, individual Cepheid distances 
and radii can be determined with an accuracy of 
the order of 5 % if the adopted /C-band surface 
brightness-color relation is correct. 

A first calibration of this relation coming di- 
rectly from interferometrically determined angular 
diameters of Cepheid variables was presented by 
Nordgrenet al. (2002) (herafter N02). They found 
a satisfactory agreement with the FG97 calibra- 
tion, within the combined 1 a uncertainties of both 
surface brightness-color calibrations. Considering 
more closely the results from N02, an even better 
agreement is found between the Fy {V — K) rela- 
tions before the zero point if forced between the 
different colors. Before this operation, N02 found 
the relation: 

iVo = 3.956 ±0.011 -0.134±0.005(F-iV')o (1) 

that translates, after forcing the zero point to the 
average of the three selected colors, to the relation: 

iVo = 3.941±0.004-0.125±0.005(y-ii')o. (2) 

On the other hand, FG97 obtain: 

= 3.947 ±0.003 -0.131 ±0.003(y-Js:)o (3) 

From this comparison, it appears that the slope 
initially determined by N02 for Fv(V — K) is sig- 
nificantly different both from their final value and 
from the FG97 relation. This difference could 
cause a bias due to the averaging of the multi- 
color zero points. Though small in absolute value, 
such a bias is of particular importance for £ Car, 
due to its relatively large {V — K) color. 

Another argument in favor of the FG97 surface 
brightness relation is that it relies on a sample of 



11 Ccpheids with periods of 4 to 39 days, while the 
relations established by N02 were derived from the 
observations of only 3 Cepheids with periods of 5 
to 10 days. Such short period Cepheids are sig- 
nificantly hotter than ^ Car (P = 35.5 days), and 
a local difference of the slope of the IRSB rela- 
tions cannot be excluded. For these two reasons, 
we choose to retain the FG97 calibration for our 
analysis of the Cepheid £ Car in the following sec- 
tion. 

4. Diameter and distance 
4.1. Angular diameter 

We have combined the photometric data from 
Pel (1976) and Bersier (2002) to construct the V- 
band light curve for t Car. The two data sets are 
spanning almost 30 years and allow an improved 
determination of the period of this variable. We 
find P = 35.54804 days. The time of maximum V 
light has been adopted from Szabados (1989) who 
give a value of Tq = 2440736.230 which is also in 
good agreement with the more recent data. The 
resulting light curve is shown in Fig.l. The K 
band light curve is based on the data from Laney 
& Stobie (1992) and is also shown in Fig.l. The 
{V — K) color curve which is needed by the IRSB 
method has been constructed on the basis of the 
observed V band data and a Fourier fit to the K- 
band data as described in Storm et al. (2004). 

For the radial-velocity curve we have used 
the data from Taylor et al. (1997) and Bersier 
(2002). Using the new ephemeris from above we 
detected a slight offset of 1.5 km s~^ between the 
two datasets. We choose to shift the Taylor et 
al. (1997) dataset by —1.5 kms^^to bring all the 
data on the well established CORAVEL system 
of Bersier (2002). We note that the exact radial 
velocity zero point is irrelevant as the method 
makes use of relative velocities. The combined 
radial velocity data are displayed in Fig.l. 

The application of the IRSB method has fol- 
lowed the procedure described in Storm et al. 
(2004). We have adopted the same reddening law 
with Rv = 3.30 and Rk = 0.30, a reddening of 
E{B -V) = 0.17 (Fernie 1990), and a projec- 
tion factor, p, from radial to pulsational velocity 
of j3 = 1.39 - 0.03 log P = 1.343 (Hindsley & Bell. 
1986; Gieren et al. 1993). As discussed by Storm 
et al. (2004) we only consider the points in the 



3 



phase interval from 0.0 to 0.8 (phase zero is de- 
fined by the V band maximum fight). We have 
appHed a smafi phase shift of —0.025 to the radial 
velocity data to bring the photometric and radial 
velocity based angular diameters into agreement. 
We note that a similar phase shift can be achieved 
by lowering the systemic velocity by 1.5 kms""*^ . 

The angular diameter curve obtained from the 
photometry has been plotted in Fig. 2, together 
with the linear displacement curve. The photo- 
metric and interferometric diameter curves are di- 
rectly compared in Fig. 3, where they are plotted 
as a function of phase. With these data we can 
compute the average angular diameters obtained 
from each technique. For the IRSB, we find an 
average limb darkened angular diameter ^?ld = 
2.974 ±0.046 milliarcsecond, and for the interfero- 
metric measurements we find ^ld = 2.992 ± 0.012 
milliarcsecond. The agreement between these two 
values is strikingly good. This is a serious indica- 
tion that the calibration of the surface brightness- 
color relation (FG97), based on non-pulsating gi- 
ant stars, does apply to Cepheids. 

4.2. Distance 

The surface brightness method yields a distance 
of 560±6 pc, and a mean radius of iZ = 179±2 Rq . 
The corresponding mean absolute V magnitude is 
My = — 5.57±0.02 mag and the distance modulus 
is (m — M)o = 8.74±0.05. The error estimates are 
all intrinsic la random errors. In addition to these 
random errors, a systematic error of the order of 
4% should be taken into account, as discussed by 
Gieren et al. (1997). The final IRSB values are 
thus d = 560 ± 23 pc, and i? = 179 ± 7 i?© . 
Compared to Storm et al. (2004) we find a sig- 
nificantly (0.24 mag) shorter distance modulus for 
£ Car. This can be explained by the use in the 
present Letter of the new and superior radial ve- 
locity data from Taylor et al. (1997) and Bersier 
(2002). 

K04 found d = 603tlt pc, using the interfero- 
metric angular diameters and a subset of the radial 
velocity data used here. To make the comparison 
more relevant, we determined the distance and ra- 
dius using the same data (interferometric diame- 
ters from Table 1 and radial velocity from Tay- 
lor et al. (1997) and Bersier (2002) - see above), 
the same ephemeris, and the same projection fac- 
tor (see Sect 4.1). Using the method of K04, we 



find a distance d = 566_;^g, and a linear radius 
R = 182tr Rq . This is in excellent agreement 
with the values obtained from the IRSB method. 

This 6% difference in the distances based on 
interferometric diameters (603 pc for K04 versus 
566 pc here) has two major causes. First, the p- 
factor used in the present paper is ~ 1.3% smaller 
than in K04. The choice of the reference used for 
the p-factor has currently an impact of a few per- 
cents on its value. This indicates that the average 
value of the p-factor for a given Ccphcid is cur- 
rently uncertain by at least a similar amount, and 
this systematic error translates linearly to the dis- 
tance determination. 

Secondly, the use of a different - and superior 
- data set for the radial velocity makes the radius 
curve different from K04. In particular the ampli- 
tude is smaller here than in K04 by ~ 3%. This 
is likely due to the more complete phase coverage 
that we have here, and possibly also to a different 
choice of spectral lines to estimate the radial veloc- 
ity. This amplitude difference translates linearly 
on the distance through the BW method. 

5. Conclusion 

The main point of our paper is to show that 
with a consistent treatment of the data, the inter- 
nal accuracy of both methods (IRSB or interfer- 
ometry) is extremely good: the angular diameter 
variation observed using the VLTI agrees very well 
with that derived from the Fv{V — K) version of 
the IRSB technique as caUbrated by FG97. For all 
the interferometric measurements, the correspond- 
ing IRSB angular diameter at the same phase lies 
within the combined Icr error bars of the two mea- 
surements (Fig. 3). Even more importantly, the 
mean angular diameter of the Cepheid as derived 
from both independent sets of angular diameter 
determination are in excellent agreement, within 
a few pcrccnts. 

Unfortunately, this is not equivalent to say that 
the Cepheid distance scale is calibrated to a 1% 
accuracy. We have drawn attention to remain- 
ing sources of systematic errors that can affect 
Cepheid radii and distances up to several percents. 
As an illustration of these sources, K04 obtain a 
distance d — 603l^g pc for £ Car, while we obtain 
d = 566^19 pc from the same interferometric data. 

We have already shown that most of the 6 % dif- 
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ferencc (equivalent to 1.3 cr) can be explained by 
the use of different radial velocity data and projec- 
tion factor. Another thing to consider is the phase 
interval used. K04 used measurements over the 
whole pulsation cycle whereas in the IRSB tech- 
nique, one avoids the phase interval 0.8-1 (Fig. 2). 
During that phase interval, that corresponds to 
the rebound of the atmosphere around the min- 
imum radius, energetic shock waves are created. 
As discussed by Sabbey ct al. (1995), they pro- 
duce asymmetric line profiles in the Cepheid spec- 
trum. Recent modeling using a self-consistent dy- 
namical approach also show that the t = 1 photo- 
sphere may not be comoving with the atmosphere 
of the Cepheid during its pulsation, at the 1 % level 
(Nardctto ct al. in preparation). Such an effect 
would impact the p-factor, modify the shape of 
the radial velocity curve, and thus bias the ampli- 
tude of the radius variation, possibly up to a level 
of a few percents. As the BW method (either its 
classical or interferometric versions) relies linearly 
on this amplitude, a bias at this level presently 
cannot be excluded. 

The interferometric BW method is currently 
limited to distances of 1-2 kpc, due to the limited 
length of the available baselines. The IRSB tech- 
nique on the other hand can reach extra-galactic 
Cepheids as already demonstrated by Gicrcn et 
al. (2000, for the LMC) and by Storm et al. (2004, 
for the SMC). Using high precision interferometric 
measurements of £ Car and other Cepheids, it will 
be possible to calibrate the IRSB method down 
to the level of a few percents. From the present 
comparison, we already see that this fundamental 
calibration will be very similar to the calibration 
found by FG97 and Nordgren et al. (2002). 

WPG acknowledges support for this work from 
the Chilean FONDAP Center for Astrophysics 
15010003. 
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Table 1: Angular diameter measurements of t Car. 

The statistical and systematic calibration uncer- 
tainties arc mentioned separately in brackets. 
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This 2-column preprint was prepared with the AAS lAlJ^X 
macros v5.2. 
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Fig. 1. — Radial velocity curve of (. Car (upper 
panel) using data from Taylor et al. (1997) shifted 
by —1.5 kms~^ (circles) and from Bersier (2002) 
(triangles). The K band photometric measure- 
ments (middle panel) were taken from Lancy & 
Stobie (1992). We have rehed on Pel (1976) (cir- 
cles) and Bersier (2002) (triangles) for the V band 
data. 
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Fig. 3. — The interferometrically determined 
angular diameters, plotted against phase (solid 
dots) with the angular diameters derived with 
the infrared surface brightness method overplot- 
ted (crosses). In the upper right corner a typical 
error bar for the surface brightness method data 
is shown. 
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